Squeezing and robustness of frictionless cooling strategies 
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Quantum control strategies that provide shortcuts to adiabaticity are increasingly considered in 
various contexts including atomic cooling. Recent studies have emphasized practical issues in order 
to reduce the gap between the idealized models and actual ongoing implementations. We rephrase 
here the cooling features in terms of a peculiar squeezing effect, and use it to parametrize the 
robustness of frictionless cooling techniques with respect to noise-induced deviations from the ideal 
time-dependent trajectory for the trapping frequency. We finally discuss qualitative issues for the 
experimental implementation of this scheme using bichromatic optical traps and lattices, which seem 
especially suitable for cooling Fermi-Bose mixtures and for investigating equilibration of negative 
temperature states, respectively. 
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I. INTRODUCTION 

Quantum control techniques that provide shortcuts 
to adiabaticity, recently proposed in the atomic physics 
framework [l[, are experiencing rapid growth due to the 
necessity, in various contexts, to reduce the undesired 
effects of a prolonged exposure to the external environ- 
ment. Such technique produces the same net result as if 
a full adiabatic following has been implemented, but on 
a much shorter timescale. More specifically, in a short- 
cut to adiabaticity technique the composition of eigen- 
states one started with is preserved at the end of the 
process, but not during the process, unlike full adiabatic 
following. The first experimental implementations have 
demonstrated efficient fast decompression of ^^Rb atoms 
in normal [2] and Bose-condensed 3 states, as well as 
fast atomic transport Q . Consequently, a number of con- 
tributions have focused attention on possible deviations 
from ideal realizations which could constitute bottlenecks 
to further implementations of these techniques. In one 
study, it was found that the transient energy excitation 
limits the maximum attainable speed [sj. In another, de- 
viations from harmonicity in the trapping potential, with 
specific attention to the crosstalk between the dynamics 
in different trapping directions, have recently been dis- 
cussed [1]. 

Our previous work showed that frictionless cooling 
(also known as "fast adiabatic" cooling), an example of 
shortcut to adiabaticity technique, may be optimal in 
controlling the temperature of a buffer gas in the sym- 
pathetic cooling of another species, such as a Fermi gas. 
This is because this cooling method allows for retaining 
the maximum value of the gas heat capacity, both as the 
gas does not enter the degenerate regime, and as it avoids 
the intentional loss of atoms required in evaporative cool- 
ing. Even in this case, we found that practical issues may 
limit the successful application of this technique. Specif- 
ically, it was shown that, for short cooling times, atoms 



temporarily attain temperatures higher than the initial 
temperature before settling down to the final tempera- 
ture. This can potentially result in atom losses in any 
realistic potential with a finite trapping depth. Also, 
very small final temperatures may result in a large spa- 
tial mismatch between the sizes of the two clouds making 
sympathetic cooling less effective. These considerations 
put practical restrictions on this method. 

Here, we extend our analysis by studying the robust- 
ness of the frictionless cooling method with respect to the 
likely presence of noise in a potential experiment. This 
is in light of the fact that, experimentally, there may be 
systematic sources of error which are of a truly nonlinear 
nature, such as noise in the trapping frequency or anhar- 
monicities in the trapping potential. To model the noise, 
we assume that an experimentalist has a well-defined pro- 
tocol for implementing frictionless cooling and that the 
effect of various noise sources present in the apparatus 
(such as current noise in the coils for magnetic trapping, 
or laser power fluctuations in an optical dipole trap) add 
up to deviations from the ideal (Ermakov) trajectory. 

In order to parametrize the robustness, one needs to 
estimate how much the eigenstate composition of the fi- 
nal state deviates from that of the initial state. However, 
usual measures such as the fidelity of the wave function 
F — \{'ip{ti)\'4t{Q))\ are unsuitable in this scheme since, 
due to the trap relaxation, the final wave function spreads 
out in space significantly. To circumvent this constraint, 
we note that our particular frictionless cooling methodQ 
can be formally linked to squeezing, and this allows us 
to use the amount of squeezing produced as a measure 
of robustness. In particular, starting from a minimum 
uncertainty state, frictionless cooling should result in a 
minimum uncertainty state at the end of the run. Any 
deviation may be considered an indicator of less than 
perfect recovery of the desired final state, i.e. a sign of 
an imperfect "adiabatic following." It should be noted 
that there are limitations in using squeezing as a gen- 
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eral measure of robustness in situations other than the 
specific coohng scheme under consideration. Specifically, 
one cannot easily generalize this measure to more so- 
phisticated states, such as mixed states or non-minimal 
uncertainty states. As a supplementary figure of merit of 
robustness, we have also chosen a modified fidelity mea- 
sure T' = \{i}(tf)\\litarget)\ whcrc and ■\l)target are 
the final states attained in the presence and in the ab- 
sence of noise, respectively, and present results using this 
measure along with the squeezing parameter. 

The paper is organized as follows: In Section II we 
discuss the connection of frictionless cooling to squeezing 
which we use as a useful indicator of the fidelity in reach- 
ing a targeted final temperature. In Section III we apply 
these considerations to discuss the robustness of the ideal 
shortcut to adiabaticity trajectories to noise fluctuations 
in the trapping parameters. The possible noise sources 
are characterized as broadband and monochromatic, and 
we discuss the interplay between their spectra in relation 
to the intrinsic spectral content of the frictionless cooling 
protocol. We flnally provide various connections to recent 
developments in the control of ultracold atoms, including 
a qualitative discussion of the realization of frictionless 
cooling in a bichromatic optical trap, a promising can- 
didate for implementing the antitrapping stage required 
for very fast cooling, provided that noise fluctuations in 
the trap are kept under control. 



II. SQUEEZING IN FRICTIONLESS COOLING 

As discussed in Ref. frictionless cooling relies on 
specific dynamics of the harmonic trapping frequency to 
transfer a quantum state such as an equilibrium thermal 
state in a trap of initial trapping angular frequency 
to that of final angular frequency wt ^ wq- This en- 
sures that the energy level spacings in the final trap are 
much smaller than those of the initial trap, hence low- 
ering the overall temperature. As we have discussed in 
detail in our previous work Q , we have additionally used 
the momentum variance as an indicator of temperature 
and have shown that the final wave function has its mo- 
mentum variance reduced by a factor of /wf compared 
to that of the initial wave function. The use of the mo- 
mentum variance as an indicator of the temperature in 
an ultracold atomic cloud has been recently validated in 
dedicated experiments [I]- 

The shortcut to adiabaticity is accomplished using a 
method originally studied by Lewis and Riesenfeld [9l.llG|. 
who introduced the invariant operator related to the har- 
monic oscillator with an arbitrary time-dependent har- 
monic potential, /(t) = (6p — m&g)^/2m -I- mwQg^/(26^). 
A key feature is that the eigenstates of can be made 
to coincide with the eigenstates of the harmonic oscillator 
Hamiltonian, particularly at the beginning and at the end 
of the cooling process. Here is a time-dependent fre- 
quency scaling factor which satisfies the Ermakov equa- 
tion 6(t) 4-a;(t)^& = uj^jh^ . This can be solved by impos- 



ing boundary conditions on and its first and second 
time derivatives, while assuming both a targeted final 
trapping frequency Wf , and a total time duration for the 
protocol, if This results in a well-defined Ermakov 
trajectory for the trap frequency By constructing a 
harmonic oscillator with the Ermakov trajectory for its 
time-dependent trapping frequency, a ground state expe- 
riences dynamic changes in the trapping frequency that 
eventually results in its cooling within the specified time 
if. In this paper, we shall denote the Ermakov trajec- 
tory as in order to distinguish it from a general 
time-dependent trap frequency. 

Formally, the above method may be understood as a 
transformation of a harmonic oscillator Hamiltonian to 
the space of the invariant operator /(t). In particular, 
the net effect can be considered as an added nonlinear- 
ity to the Hamiltonian that produces squeezing of atomic 
states. A first step in this direction has been reported in 
[Tl| where number (Fock) states during sudden and adi- 
abatic expansions of the trap have been discussed. One 
can show how squeezing is a natural byproduct of the Er- 
makov dynamics by noting the relationship between the 
raising and lowering operators of the harmonic oscillator 
Hamiltonian a^(fy and a(t\ and the Dirac raising and 
lowering operators for the Lewis and Riesenfeld invari- 
ant operator /(i) a^if) and a(t)[l0|: 

= 77(t)a(t) + CWa^(t) (1) 
a\t) = CW*a(t)+^W*a^W (2) 



where 




, (3) 



(4) 



Despite the complicated expressions for 77 (t) and C(Oi 
Eqs. (m and ^ are essentially Bogoliubov transforma- 
tions, well-known in the description of squeezing in quan- 
tum optics [T^]. In particular, and C(0 satisfy the 
condition |?7p — |CP = li one can directly recover the 
standard squeezed state result of quantum optics [l2| and 
write a = 5a5^ where 5 is a unitary squeezing operator 
S{i) = exp(-^at2/2-heaV2) with ^ = cosh'^ \T]{t)\ 
which, acting on a vacuum, produces a squeezed state 
10 = S'(O|0). This connection went unnoticed in the 
paper by Lewis and Riesenfeld [loj . presumably because 
the field of quantum optics was still in its infancy. 

While the position variance of the wave function un- 
dergoing the frictionless cooling is analytically given as a 
function of the solution to the Ermakov equation b{t) [l| 

it can be shown that the momentum variance, propor- 
tional to the temperature of the gas, evolves as 

hmuj[ 



■exp(2cosh"^ \r]{t)\), 



(6) 



so that their product is 

4itHit) = ^^b'it)exp{2cosh~' by(t)l). (7) 
4 Wo 

This shows that the frictionless coohng strategy formally 
produces a squeezed state, and may be potentially useful 
as a quantum control method for generating s queezin g in 
cold atoms. At the final time t = , b{tf) = \J wq /wt and 
exp(2cosh"^ l?7(^f)l) = 1, such that a^(tf)cr2(tf) = W-j^ 
i.e. the state returns to the minimum uncertainty state as 
desired. This implies that one may quantify the robust- 
ness to external noise by checking how much the state 
differs from the minimum uncertainty state when the 
time-dependent trap frequency deviates from the opti- 
mal Ermakov trajectory. 

To investigate quantitatively the squeezing effect in 
frictionless cooling strategies, we have numerically solved 
the Schrodinger equation for the time-dependent har- 
monic oscillator with the trap frequency given by ujE{t). 
As a starting point, we have compared the two decom- 
pression strategies for a harmonically trapped atomic 
cloud we have used in our previous work [7[ where, for 
concreteness, we assumed the initial trap angular fre- 
quency Wo = 27r X 250 Hz and the final trap angular fre- 
quency Wf = 27r X 2.5 Hz. The two strategies are charac- 
terized by their duration: if =25 ms for which uj%{t) > 
throughout the entire run, and tf=6 ms which includes a 
short time interval for which uj%{t) < 0. During this in- 
terval the atoms temporarily experience an antitrapping 
potential (like in an inverted harmonic oscillator, with 
the possibility of cooling first noticed in the conclusions 
of [iJl) leading to an accelerated spreading of the wave 
function, after which the trap is flipped back to allow 
regrouping, and desired cooling, of the wave function by 
t = t{. 

We show in Fig. 1 the time dependence of the position 
and momentum variances as well as their product for 
the two durations of cooling. The temporal variation of 
the position variance is found to be independent of if, 
as expected from the analytical results, while the time 
evolution of the product of variances clearly indicates 
the recovery of the minimum uncertainty state at t = t{. 
For comparison, a quasi-adiabatic trajectory obtained by 
a linear ramping down of the frequency to its final value 
is also depicted, showing that the uncertainty product in 
this case exceeds the minimum value at final time. 

It should be noted here that, since Eqs. ([5])- ([7]) are 
explicitly given in terms of the Ermakov solution b{t), it 
is possible to directly calculate the amount of squeezing 
without evolving the wave function using the Schrodinger 
equation with a time-dependent potential. Figure 2 
shows the momentum variance and the product of po- 
sition and momentum variances as a function of time for 
various possible values of tf. The results for t{=6 ms and 
if =25 ms obtained by evolving the Schrodinger equation, 
were found accurately reproduced. It is clear that the 
methods of quantum optical squeezing are applicable in 
the context of frictionless cooling, although the system is 
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FIG. 1: (Color online) Squeezing properties of the frictionless 
cooling strategies. Depicted is the time dependence of the po- 
sition variance (top), momentum variance (center), and their 
product (bottom) for the three cooling strategies discussed 
in the text, with the product cr^cTp expressed, from now on, 
in units of . The dashed and solid lines are for frictionless 
cooling with tf =6 ms and 25 ms, respectively. The dot-dashed 
line is for a linear ramp-down of the frequency corresponding 
to a quasi-adiabatic expansion occurring in a time ti = 2Ti/uJt 
as discussed in 01 • The adiabatic invariant strategies allow 
for the preservation of states with minimal uncertainty prod- 
uct, as they display an uncertainty product returning to the 
value (j^(Tp = at t = t{ after a transient increase at ear- 
lier times. Conversely, the quasi-adiabatic expansion keeps 
its uncertainty product constant for a long time interval, but 
manifests an increase approaching the final time as the quasi- 
adiabaticity condition is less satisfied in the latest stage with 
minimal trapping frequency. 

rather different from the standard quantum optical sys- 
tems involving nonlinear optical media. 

For small if, it is noted that the momentum variance 
<Jp{t) goes up significantly during the run, resulting in 
large values of the variance product, a'^{t)ap{t). The 
presence of such large deviation from the minimum un- 
certainty state during the run implies that, for such cases, 
any noise is likely to more effectively push the system 
away from the minimum uncertainty state. Indeed, we 
observe such a trend in the next section. 



III. ROBUSTNESS OF THE FRICTIONLESS 
COOLING SCHEME 

In this section, we want to answer the following ques- 
tion: how robust is the frictionless cooling technique to 
small deviations from the optimized Ermakov trajectory? 
Since, as discussed in the last section, the variable fre- 
quency strategy is basically mappable to a nonlinear sys- 
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FIG. 2: Top: logj^,, calculated directly from the squeez- 

ing transformation for various t[. Bottom: logj^o[a^(f)(Tp(t)] 
for various t{ indicating the deviation from, and the eventual 
return to, the minimum uncertainty state over time. It is 
noted that for small if, the momentum variance (Jp{t) pro- 
portional to the temperature, goes up significantly during the 
run. As previously noticed in our work Q, due to potential 
atom loss, this presents a restriction on using this method in 
sympathetic cooling. A consequence relevant to this paper is 
the very high values of the product a'^{t)ap(t) during the run 
for short tf. 

tem, the interplay between the deterministic Ermakov 
trajectory and any added noise may, in principle, be sub- 
tle and can lead to large and detrimental deviations from 
the idealized behavior. To study robustness in this sys- 
tem, we consider two experimentally relevant types of 
noise, broadband and monochromatic. These choices al- 
low us to understand the system via the complementary 
responses. 

A. Response to random noise 

For the random noise, we consider two different possi- 
bilities: Gaussian white noise with mean value lo'^ (t) , and 
standard deviation /3|a;|;(t)|, and uniformly distributed 
noise centered at with width /3\Ld'^{t)\, where (3 is 
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FIG. 3: (Color online) The final variance product, 
(7^(tf)(7p(tf) (bottom solid line, left vertical axis) and the fi- 
delity (top dashed line, right vertical axis) for t{=25 ms (red 
squares) and ti—6 ms (black circles) as a function of noise 
amplitude /3 described in the text. Gaussian white noise (top 
plot); uniform white noise (bottom plot), all from the averag- 
ing of 150 realizations. 



varied from zero to 0.2, i.e. in both situations noise has 
an amplitude up to 20% of the unperturbed trajectory 
|w|;(i)|. The results of simulations with the Gaussian 
white noise and uniform white noise are shown in Fig. [3] 
where the final variance product, al{tf)<7p{tf) and the fi- 
delity, obtained by averaging 150 trajectories, is depicted. 
Both the variance product and the fidelity show that the 
overall deviation from the desired result is of order 0.1%, 
and hence the system is relatively insensitive to random 
noise. This scheme is therefore more robust than we ex- 
pected, especially given the very large amount of noise 
(up to 20%) added. This may be a result of the versatil- 
ity of the Ermakov invariant scheme which, in principle, 
is able to handle all values of tf; the net effect of added 
noise may be viewed as shifting the final t{ which still 
produces the correct state. As expected, the final state 
deviates more and more from the desired state as one 
increases the noise amplitude (3. It is worth noting that, 
as we had expected from examining Fig. [21 the if =6 
ms case with negative square frequency in the Ermakov 
trajectory is found to be more sensitive to the effect of 
noise. Physically, this is reasonable considering the added 
vulnerability to noise during the time interval in which 
atoms experience their antitrapping stage. 
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FIG. 4: Temporal density of states versus uj'^ /ujq for the two 
adiabatic strategies discussed in the text with t{=25 ms (top) 
and 6 ms (bottom). The temporal density of states, here 
evaluated for 10^ samples and 10"^ squared frequency intervals 
in the entire range of variability of uj'^ , is normalized to unity, 
thereby representing the probability density to realize a given 
Lo'^ during the corresponding frictionless cooling process. 



B. Response to sinusoidal modulation 

Noting that realistic random noise can be viewed as a 
sum over sinusoidal modulations with random frequen- 
cies and amplitudes, one can analyze the robustness of 
this system by studying the response to a range of purely 
sinusoidal modulations. Sinusoidal modulations may oc- 
cur as an unintended noise or added intentionally, for 
instance to improve the signal-to-noise ratio of a specific 
feature under investigation, to create arbitrary trapping 
potentials or to prepare well-defined target states, 
for instance atomic Fock states in optical traps [TB - flsj 
and optical lattices [l^, [2l| . Continuing our analogy to 
quantum optics, the study of monochromatic noise is 
similar to examining the response of a nonlinear optical 
medium to coherent light of fixed frequency. 

Before such analysis we need to speciiy the range of 
trapping frequencies involved in the cooling procedure. 
Wc indeed expect that the effect of noise sources will de- 
pend on their relative location with respect to the spec- 
trum of frequencies of the system. To quantify the latter. 



we introduce a 'temporal' density of states defined as the 
amount of time spent by the system in the range be- 
tween wl; and w|; + d{oj'^). This is directly related, at 
least for monotonic decreases of the trapping frequency, 
to [d{uj'^)/dt]^^ , the inverse of the slope of the curve 
uj'^(t). For the case of nonmonotonic behavior, such as 
in the strategics involving a stage in which < 0, this 
is less trivial to obtain analytically. However, in practice, 
we can use a simple code in which we sort the times at 
which frequencies occur and count the number of times 
we find a square frequency in the [ui'^^uj'^ + d{u!%)] in- 
terval. For the frictionless cooling strategy with tf=25 
ms this is depicted in the top plot of Fig. 4. The case 
of strategies in which there are also time intervals with 
negative square frequencies is more involved, and a typi- 
cal example is shown in the bottom plot of Fig. 4 for the 
t[—6 ms case. Multiple peaks now occur - one located 
at the minimum negative square frequency, as the system 
spends a large amount of time in that region, and another 
at a positive square frequency where the plot shows a lo- 
cal maximum in the final stage of cooling. The peak 
at zero frequency gets contributions from three distinct 
time intervals, two with a zero-crossing behavior, before 
and after reaching the minimum square frequency, and 
one from positive square frequencies alone at the very 
end of the cooling (sec Fig. 1 in 0), resulting in a dis- 
continuity in the temporal density of states. The high 
peaks in the < region for the t[ =6 ms case indi- 
cate that the contribution from the antitrapping stage is 
much more significant than what a simple inspection of 
the Ermakov trajectory would suggest. 

The effect of sinusoidal noise was studied by modify- 
ing the Ermakov trajectory u>'^{t) to uj'^E{t) — w|;(t)[l 4- 
/? sin(r2„t)] and evolving the Schrodinger Equation under 
this time-dependent potential. In Fig. [5]wc plot the fi- 
nal variance product, (T^(if)o'p(if) and the fidelity (inset) 
for various tf and (3. We have chosen two possible cases 
of the modulation frequency, f2„ < ojq (0„ = ujq/2) and 
fi„ > Wo (f^n = 3a;o/2). The results show that, as to be 
expected; the deviation from the desired state increases 
with increasing /3, although the effect is now much more 
pronounced than those for the random noise. Again, as 
we expected from Fig. [21 it is seen that the short tf cases 
that involve negative square frequencies are more sensi- 
tive to noise. There is also an overall oscillatory behavior 
as a function of tf for the variance product, and this may 
be explained as follows. If one assumes, as a very crude, 
first order approximation, that the only relevant quan- 
tity in the wave function evolution is the instantaneous 
trapping frequency in the presence of noise, i.e. we com- 
pletely disregard the trajectory history, one can write 
using Eq. ([7]) the final position variance 

^I(ify^(if) = f ^, (8) 

4 Wf 

where ujf is the new final trap frequency in the presence 
of the added noise. Given the form of the noise, wc have 



t, (ms) t| (ms) 

FIG. 5: The final variance product, a^{ti)ap{ti) and the fidehty (inset) as a fiinction of tf and the amplitude of sinusoidal 
modulation in a, modulated trajectory LUs{i) = + /3sin(n„t)]. In the left panel we show the case of a sinusoidal noise 

within the range of trapping frequencies, Qn ~ i^o/2, in the right panel an example of noise outside, Q„ = 3(j;o/2. 



(Df = cjf-^/l + /3 sin(ri„if ), and hence 

^littHitt) - ^Vl + /5sin(r!„tf), (9) 

so that there is an oscillatory behavior expected for the 
variance product as a function of tf at constant f2„ (and 
vice versa). This result indicates a fairly strong depen- 
dence of robustness on t{ or, to be more precise, on the 
deviation from the desired final trap frequency due to 
the presence of noise which, in this case, is given by the 
multiplicative factor + /3 sin(fi„tt). 

Next, we have fixed the amplitude of the sinusoidal 
noise (3 to be 0.1 i.e. 10% of the instantaneous Er- 
makov trajectory u!%{t), and calculated the response of 
the system to noise in a range of modulation frequencies, 
rin G [0,2a;o]. This is presented in Fig. |6] for both the 
if =25 ms and if =6 ms cases. It is observed that the sys- 
tem is fairly sensitive to the variation of the modulation 
frequency, and it is also evident that the tf= 6 ms case 
with the negative square frequency stage shows a very 
strong deviation from the desired state. As before, this 
is likely due to the enhanced vulnerability to external 
modulation during the anti-trapping stage. The oscilla- 
tory behavior as a function of the modulation frequency 
rin for fixed if may, again, be partially explained using 
Eq. ^ . From the same equation, the maximum possible 
deviation from the desired state is expected to be of the 
order 13/2. Since we use /3= 0.1 for the maximum noise, 
this rough estimate holds for the if=25 ms case but not 
for the if=6 ms case. The fact that the maximum pos- 
sible deviation in variance product worked out from Eq. 
(|9]) is smaller than that observed indicates the presence 
of other contributing factors that disrupts the shortcut 
to adiabaticity as one would expect. 



In order to produce a strong response, the modulation 
frequency n„ should, in general, match those trap fre- 
quencies along the trajectory that varies slowly i.e. those 
with large dt/d{uP) in Fig. 2] For if = 25 ms this seems 
to be the case, in that the strongest response appears near 
ri„ sa 0. However for if = 6 ms, although a strong re- 
sponse is also visible near the expected region of 17„ sa 0, 
the maximum peak clearly appears near il^/wg ~ 0.5. 
On close examination, it was found that the peaks ap- 
pear near certain multiples of a characteristic angular fre- 
quency wc = 7i'/2if, obtainable from Eq. ^ as the first 
local maximum. It was observed that the the first three 
peaks are found at O = ujc, 'i^Ci and 7ujc. Comparing 
the square of characteristic angular frequency, the if = 25 
ms case gives fl^/uj^ = 0.0016, 0.0256, 0.0784, while the 
if = 6 ms case gives f^V^o = 0.0278,0.4444, 1.3611, ap- 
proximately coinciding with the corresponding locations 
of the maximum peaks. This shows that, in addition to 
the temporal density of states, the characteristic angu- 
lar frequency that depends on the precise type of noise 
and its effect on the final trapping frequency is another 
important parameter in characterizing the robustness of 
this system. 

Finally, wc note that the overall behavior of the fidelity 
mirrors that of the variance product, indicating that, in 
our case, the amount of squeezing is indeed an efficient 
alternative paramctrization of robustness. 



C. Reduction of the maximum temperature 

A crucial finding in our previous work Q was the lim- 
itation imposed by the high temperature attained by the 
atoms during the cooling process. Therefore, related to 





FIG. 6: (Color online) The variance product, <Jx{ti)o'p(t() 
(bottom solid line in black , left vertical axis) and the fi- 
delity (top solid line in red , right vertical axis) as a function 
of the scaled modulation frequency SI^/cjq for fixed modula- 
tion amplitude P — 0.1 with tf=25 ms (top plot) and tf=6 ms 
(bottom plot). 




FIG. 7: Scaled maximum temperature as a function of the 
the scaled modulation frequency (Q'^/luq) for fixed modula- 
tion amplitude /3 — 10% with tt—25 ms (top) and t{=6 ms 
(bottom). 



the study of robustness of the system to noise, we have 
investigated the maximum temperature attained by the 
atoms during the cooling process in the presence of si- 
nusoidal modulation of various modulation frequencies. 
For the Ermakov trajectory which does not contain an 
antitrapping stage, the sinusoidal modulation can pro- 
vide momentum kicks, increasing the momentum vari- 
ance. On the other hand, for the t{—6 ms case that in- 
cludes an antitrapping stage, with a modulation of the 
curvature of the anti-trap one may be able to "catch" 
the wave function spread in such a way that the mo- 
mentum variance is reduced. Qualitatively, this may be 
viewed similarly to the mechanism behind the dynamic 
equilibrium attained by an inverted pendulum with an 
oscillating pivot. Figure [7] shows the scaled maximum 
temperature as a function of the modulation frequency. 
It is found that for the ti~25 ms case, there is some heat- 
ing as the modulation frequency is increased. However, 
interestingly, for the tf =6 ms case that contains the neg- 
ative square frequency stage, it shows a dip in the max- 
imum temperature around the modulation frequency of 
rin = l.SSwo to a value less than that corresponding to 
no sinusoidal modulation (fi„ = 0). In addition, the vari- 
ance product at this frequency is also very close to the 



minimum uncertainty state, indicating this is the opti- 
mum set-up for the tf= 6 ms case if reducing atom losses 
is the most important priority. This provides a motiva- 
tion for intentionally adding a sinusoidal modulation in 
the cooling scheme to minimize atom-loss due to heating, 
adding a further dimension to the issue of cooling opti- 
mization when an antitrapping stage is present [22| . It 
has been shown [23| that the cooling time is minimized, 
in the presence of a trapping stage, using a "bang-bang" 
control, i.e. with stepwise changes of the trapping fre- 
quency, but the possible addition of a sinusoidal driving 
may call for a modification, or at least a check, of this 
optimization procedure. 



IV. DISCUSSION 

The results reported in this paper are encouraging for 
the experimental implementation of frictionless cooling, 
possibly including a negative square frequency stage, in 
Fermi-Bose mixtures using optical dipole traps and op- 
tical lattices. A scheme using a single-frequency optical 
dipole trap with a continuously tunable laser is not feasi- 
ble due to the large atom losses expected in crossing the 
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dominant atomic transition from the red-detuned to the 
blue-detuned side to achieve antitrapping. Nevertheless, 
this issue may be circumvented through the use of two 
laser beams at constant frequencies, opposite detunings, 
and variable power ratio, such as the bichromatic optical 
dipole traps proposed in Q to allow optimal heat capac- 
ity matching between the Bose and Fermi gases [25l. |26{. 
The presence of noise related to power fluctuations and 
beam-pointing stability for both laser beams in this con- 
figuration makes the discussion presented in this paper 
quite relevant for implementing frictionless cooling in 
bichromatic optical dipole traps. Bichromatic trapping 
has been recently implemented at the magneto-optical 
trap stage for a single species [l^], and two-species se- 
lective trapping and cooling has been demonstrated with 
hybrid traps involving magnetic and optical confinement 
[28I [29} . Therefore experiments involving trapping and 
cooling of two species in a bichromatic trap should be 
within reach in the near future. Our discussion should 
also be relevant to the case of frictionless cooling in opti- 
cal lattices, as recently discussed [13] ■ Dynamically vari- 
able spacing (the so-called accordion lattice) allows for a 
continuous increase of lattice periodicity, as experimen- 
tally demonstrated in (3ll - [33| , without the need to change 
the laser frequency. In this case an additional source of 
noise during the trap expansion is due to the presence of 



acousto-optic deflectors, adding up to the beam-pointing 
stability of the lasers. Implementing frictionless cooling 
with a negative square frequency stage via a bichromatic 
optical lattice could also be of great relevance to inves- 
tigate fundamental issues of statistical mechanics as the 
approach to equilibrium in atomic systems at negative 
temperatures |34l - l36l |. 

In summary, we have found that the robustness of the 
frictionless cooling method to noise can be analyzed by 
characterizing the flnal quantum state in terms of the 
amount of squeezing as well as the usual fidelity. It was 
found that the method is quite robust to the presence of 
broadband noise in the trapping frequency, and further 
analysis involving monochromatic sinusoidal modulation 
has allowed us to resolve its response. We have studied 
the dependence of the squeezing and fidelity on the de- 
viation from the expected Ermakov trajectory, and iden- 
tified the role of such measures as the temporal density 
of states and characteristic angular frequency. Numeri- 
cal simulations indicate that, despite the perception that 
short ti could mean less time for noise to interfere with 
the system, too short a t{ is best avoided in practice. In 
addition, a way to reduce the maximum temperature and 
hence atom losses was found by adding a high frequency 
sinusoidal modulation, which helps to mitigate one of the 
limitations of this scheme found in our previous work. 



[1] X. Chen, A. Ruschhaupt, S. Schmidt, A. del Campo, D. 
Guery-Odelin, and J. G. Muga, Phys. Rev. Lett. 104, 
063002 (2010). 

[2] J.-F. Schaff, X.-L. Song, R Vignolo, and G. Labeyrie, 

Phys. Rev. A 82, 033430 (2010). 
[3] J.-F. Schaff, X.-L. Song, P. Capuzzi, P. Vignolo, and G. 

Labeyrie, EPL 93, 23001 (2011). 
[4] E. Torrontegui, S. Ibanez, X. Chen, A. Ruschhaupt, D. 

Guery-Odelin, and J. G. Muga, Phys. Rev. A 83, 013415 

(2011). 

[5] X. Chen and J. G. Muga, Phys. Rev. A 82, 053403 
(2010). 

[6] E. Torrontegui, X. Chen, M. Modugno, A. Ruschhaupt, 
D. Guery-Odehn, and J. G. Muga, Phys. Rev. A 85, 
033605 (2012). 

[7] S. Choi, R. Onofrio, and B. Sundaram, Phys. Rev. A 84, 

051601(R) (2011). 
[8] M. J. Davis, P. B. Blakie, A. H. van Amerongen, N. J. 

van Druten, and K. V. Kheruntsyan, Phys. Rev. A 85, 

031604 (2012). 
[9] H. R. Lewis, Phys. Rev. Lett. 18, 510 (1967). 
[10] H. R. Lewis and W. B. Riesenfeld, J. Math. Phys. (N.Y.) 

10, 1458 (1969). 
[11] D. Sokolovski, M. Pons, A. del Campo, and J. G. Muga, 

Phys. Rev. A 83, 013402 (2011). 
[12] D. F. Walls and G. J. Milburn, Quantum Optics 

( Springer- Verlag, Berlin, 1994); S. M. Barnett and P. M. 

Radmore, Methods in Theoretical Quantum Optics (Ox- 
ford University Press, New York, 2002). 
[13] C. Yuce, A. Kilic, and A. Coruh, Phys. Scr. 74, 114 

(2006). 



[14] Ph. W. Courteille, B. Deh, J. Fortagh, A. Giinther, S. 

Kraft, C. Marzok, S. Slama, and C. Zimmermann, J. 

Phys. B 39, 1055 (2006). 
[15] A. M. Dudarev, M. G. Raizen, and Q. Niu, Phys. Rev. 

Lett. 98, 063001 (2007). 
[16] A. del Campo and J. G. Muga, Phys. Rev. A 78, 023412 

(2008) . 

[17] M. Pons, A. del Campo, J. G. Muga, and M. G. Raizen, 

Phys. Rev. A 79, 033629 (2009). 
[18] S. Wan, M. G. Raizen, and Q. Niu, J. Phys. B 42, 195506 

(2009) . 

[19] M. G. Raizen, S. P. Wan, C. Zhang, and Q. Niu, Phys. 

Rev. A 80, 030302(R) (2009). 
[20] F. Heidrich-Meisner, S. R. Manmana, M. Rigol, A. Mu- 

ramatsu, A. E. Feiguin, and E. Dagotto, Phys. Rev. A 

80, 041603(R) (2009). 
[21] G. M. Nikolopoulos and D. Petrosyan, J. Phys. B 43, 

131001 (2010). 

[22] D. Stefanatos, H. Schaettler, and J.-S. Li, SIAM J. Con- 
trol Optim. 89, 2440 (2011). 

[23] D. Stefanatos, J. Ruths, and J.-S. Li, Phys. Rev. A 82, 
063422 (2010). 

[24] R. Onofrio and C. Presilla, Phys. Rev. Lett. 89, 100401 

(2002) . 

[25] C. Presilla and R. Onofrio, Phys. Rev. Lett. 90, 030404 

(2003) . 

[26] R. Onofrio and C. Presilla, J. Stat. Phys. 115, 57 (2004). 
[27] C. Qiang, L. Xin-Yu, G. Kui-Yi, W. Xiao-Rui, C. Dong- 

Min, and W. Ru-Quan, Chin. Phys. B 21, 043203 (2012). 
[28] J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, G. 

Thalhammer, F. Minardi, S. Stringari, and M. Inguscio, 



9 



Phys. Rev. Lett. 103, 140401 (2009). 
[29] F. Baumer, F. Miinchow, A. Gorlitz, S. E. Maxwell, P. S. 
Julienne, and E. Tiesinga, Phys. Rev. A 83, 040702(R) 
(2011). 

[30] C. Yuce, Phys. Lett. A 376, 1717 (2012). 

[31] L. Fallani, C. Fort, J. E. Lye, and M. Inguscio, Optics 

Express 13, 4303 (2005). 
[32] T. C. Li, H. Kelkar, D. MedeUin, and M. G. Raizen, 

Optics Express 19, 5465 (2008). 



[33] S. Al- Assam, R. A. Wilhams, and C. J. Foot, Phys. Rev. 

A 82, 021604(R) (2010). 
[34] A. P. Mosk, Phys. Rev. Lett. 95, 040403 (2005). 
[35] A. Rapp, S. Mandt, and A. Rosch, Phys. Rev. Lett. 105, 

220405 (2010). 

[36] S. lubini, R. Franzosi, R. Livi, G.-L. Oppo, and A. Pohti, 
,arXiv:1203.4162 (19 Mar 2012). 



